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HIGHLIGHTS 


•  This  paper  presents  a  complete  and  real  structure  of  the  large-scale  PV  systems. 

•  FTA  is  used  to  analyze  the  effects  of  a  battery  system  on  the  system  reliability. 

•  We  estimate  total  component  reliability  and  overall  reliability  for  the  PV  system. 

•  Increasing  nominal  power  output  of  the  PV  system  will  decrease  Reliability. 

•  The  critical  components  with  priority  for  the  PV  system  are  revealed. 
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Photovoltaic  (PV)  systems  have  significantly  shifted  from  independent  power  generation  systems  to  a 
large-scale  grid-connected  generation  systems  in  recent  years.  The  power  output  of  PV  systems  is 
affected  by  the  reliability  of  various  components  in  the  system.  This  study  proposes  an  analytical 
approach  to  evaluate  the  reliability  of  large-scale,  grid-connected  PV  systems. 

The  fault  tree  method  with  an  exponential  probability  distribution  function  is  used  to  analyze  the 
components  of  large-scale  PV  systems.  The  system  is  considered  in  the  various  sequential  and  parallel 
fault  combinations  in  order  to  find  all  realistic  ways  in  which  the  top  or  undesired  events  can  occur. 
Additionally,  it  can  identify  areas  that  the  planned  maintenance  should  focus  on.  By  monitoring  the 
critical  components  of  a  PV  system,  it  is  possible  not  only  to  improve  the  reliability  of  the  system,  but  also 
to  optimize  the  maintenance  costs.  The  latter  is  achieved  by  informing  the  operators  about  the  system 
component’s  status.  This  approach  can  be  used  to  ensure  secure  operation  of  the  system  by  its  flexibility 
in  monitoring  system  applications.  The  implementation  demonstrates  that  the  proposed  method  is 
effective  and  efficient  and  can  conveniently  incorporate  more  system  maintenance  plans  and  diagnostic 
strategies. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Analyzing  the  reliability  of  PV  power  systems  is  important  for 
planning  and  long-term  operation,  because  the  analysis  helps 
predict  system  behavior  over  time  and  devise  appropriately  timed 
maintenance  plans.  It  is  a  significant  factor  for  the  operator  to  be 
able  to  assess  system  reliability  under  long-term  operations  in 
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order  to  optimize  decisions  in  design,  engineering,  procurement, 
construction,  and  service  [1], 

PV  and  wind  systems  produce  electric  power,  which  involves 
zero  greenhouse  gas  emissions  and  fossil  fuel  consumption.  The 
total  capacity  of  grid-connected  PV  power  systems  has  grown 
exponentially,  from  300  MW  in  2000  to  approximately  67  GW  in 
2011  [2j.  PV  power  generation  will  become  the  main  focus  of  future 
energy  development.  As  a  clean  energy,  its  application  has  gradu¬ 
ally  changed  toward  large-scale,  grid-connected  systems.  It  signif¬ 
icantly  influences  the  reliability,  economics  and  operational 
stability  of  such  systems  [3],  The  largest  PV  system  with  a  gener¬ 
ation  capacity  of  80  MW  was  installed  in  Sarnia,  Ontario,  Canada,  in 
2010  [4],  Additionally,  the  European  PV  Technology  Platform  Group 
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Nomenclature 

h 

reverse  saturation  current 

Ns 

cells  connected  in  series 

Jph 

internal  equivalent  current  source 

k 

the  Boltzmann  constant  of  1.38  x  10  23  J  K  1 

vt 

thermal  voltage  of  the  array 

a 

the  permittivity  of  the  diode 

Q 

the  electron  charge  of  1.602  x  10~19  Coulomb 

i<i 

temperature  coefficient  (mA  °C_1) 

T> 

the  temperature  of  the  pn  junction 

Cn 

the  reference  irradiation 

G 

solar  irradiation 

Eg 

bandgap  energy  of  the  semiconductor 

AT 

temperature  deviation  from  the  reference  value 

Pmax 

maximum  output  power  of  a  PV  cell 

Tn 

reference  temperature 

Vmpp.min  minimum  MPP  voltage  of  the  inverter 

Isr 

the  reference  reverse  saturation  current 

^mpp 

PV  voltage  at  the  maximum  power  point 

MPP 

maximum  power  point 

Voc 

PV  open-circuit  voltage 

fsc 

PV  short-circuit  current 

Mv 

PV  temperature  coefficient  of  Voc 

Wnpp.max 

maximum  MPP  voltage  of  the  inverter 

Wnax 

maximum  operating  voltage  of  the  inverter 

A  i 

PV  temperature  coefficient  of  fsc 

m 

reliability  probability  function 

^mpp 

PV  current  at  the  maximum  power  point 

m 

probability  density  function 

^DC.max 

maximum  input  DC  current  of  the  inverter 

GP 

grid  protection 

F(t) 

failure  probability  function 

SPD 

surge  protection  device 

ACS,  DCS 

AC  or  DC  switch 

BD 

blocking  diode 

CB 

AC  circuit  breaker 

PV 

photovoltaic  cell 

INV 

inverter 

BS 

battery  system 

CON 

connector 

DCB 

differential  circuit  breaker 

CC 

charge  controller 

1  -  Pr(£)  reliability  probability 

FV 

Fussel— Vesely 

MCS; 

minimal  cut  set  i 

Pr(£) 

failure  probability 

m 

the  number  of  basic  events  in  the  largest  minimal  cut 

Prop 

the  probability  of  top  event 

set 

Pmcs, 

the  probability  of  occurrence  of  MCS; 

n 

the  number  of  minimal  cut  sets 

Pb  1.  Pb2,- 

.Psn  the  failure  probabilities  of  basic  events 

N 

total  number  of  batteries 

M 

number  of  batteries  in  parallel 

£; 

the  event  that  component  i  operates  without  failure 

T 

the  system’s  time  to  failure 

^Charge- 

-controller  failure  rate  of  the  charge  controller 

^Battery 

failure  rate  of  the  battery  system 

reports  that  PV  systems  are  predicted  to  reach  network  parity  in 
most  of  Europe  in  2019  [5], 

In  recent  literature,  evaluating  the  reliability  of  solar  PV  has 
been  a  point  of  interest. 

In  Ref.  [6],  the  investigators  analyzed  the  reliability  of  solar  PV 
power  system  designs  using  failure  mode  effect  analysis  (FMEA) 
and  fault  tree  analysis  (FTA),  and  also  calculated  the  failure  rates  of 
the  PV  array  and  inverter.  In  Ref.  [7],  the  investigators  estimated 
reliability  equations  from  the  FTA  but  did  not  analyze  the  reliability 
probability  functions.  The  maximum  reliability  of  PV  arrays 
with  optimal  interconnection  of  PV  modules  was  investigated  in 
Ref.  [8—10],  In  Ref.  11  ],  the  researchers  evaluated  the  reliability  of 
an  electric  power  generation  system,  including  a  PV  system,  by 
considering  the  load  under  the  assumption  that  none  of  the  system 
components  ever  failed.  The  researchers  in  Ref.  [12]  proposed  a 
new  method  for  the  calculation  of  the  optimal  configuration  of 
large-scale  PV  systems.  In  Ref.  [13],  FTA  and  Markov  chain  method 
are  jointly  used  to  evaluate  the  behavior  of  PV  system.  The  energy 
cost  of  PV  system  is  estimated  and  applied  to  PV  system  designs. 
The  investigators  in  Ref.  [14]  studied  the  reliability  of  battery 
voltage  regulators  (BVRs)  used  in  PV  systems  and  calculated  the 
overall  system  reliability.  In  Ref.  [15  ,  the  investigators  proposed  a 
model  using  Monte  Carlo  for  the  analysis  of  reliability  of 
rechargeable  batteries  in  photovoltaic  power  supply  systems.  In 
Ref.  [16],  the  researchers  discussed  the  PV  inverters  used  in  PV 
systems,  presenting  their  experimental  results.  The  reliability  of  PV 
systems  was  estimated  in  various  small-scale  field  tests  described 
in  Refs.  [17—19], 

Although  a  wide  variety  of  studies  have  been  conducted  about 
the  large-scale,  grid-connected  PV  systems  [20—27],  the  real 


electrical  architecture  of  modern  large-scale,  grid-connected  PV 
systems  with  battery  backup  requires  further  consideration. 

This  paper  presents  a  technique  for  analyzing  the  reliability  of 
large-scale,  grid-connected  PV  systems  using  an  exponential  dis¬ 
tribution  based  on  FTA  method  and  considering  the  presence  of  a 
battery  system  and  charge  controller.  It  is  necessary  to  point  out 
that  changing  the  function  from  exponential  distribution  to,  for 
example,  accelerated  life  tests  (ALTs)  with  Log-normal,  Weibull  or 
mixed-WeibuIl  distributions  did  not  alter  the  output  of  the  pro¬ 
posed  method.  In  addition,  if  the  repair  interval  of  the  system 
components  is  sufficiently  less  than  a  critical  value  and  does  not 
influence  the  system  operation,  then  the  repair  time  could  be 
ignored.  Thus,  in  this  study,  it  was  assumed  that  any  failures  could 
not  be  repaired  [7],  [11],  and  [28—32];  therefore,  if  a  component 
failed,  the  overall  PV  system  was  assumed  to  be  in  failure  mode  and 
the  repair  is  finished  during  a  short  time.  However,  the  overall 
system  remaining  in  operation  during  that  time,  because  the 
electric  power  of  distribution  system  is  sufficient  for  loads  to  use.  It 
is  also  important  to  point  out  that  in  order  to  define  the  critical 
components,  the  scope  of  this  paper  is  focused  on  the  reliability 
evaluation  of  PV  systems  based  on  unrepairable  components.  It  is 
assumed  the  repair  time  of  the  system  failure  is  too  short  and  the 
loads  could  be  incorporated  with  the  distribution  system  during 
failure.  However,  the  proposed  method  is  applicable  when  the 
repair  and  the  common-cause  degradation  are  considered. 

The  remainder  of  the  paper  is  organized  as  follows.  Section  2 
describes  the  electrical  structure  of  large-scale  grid  connected  PV 
systems.  Section  3  proposes  the  reliability  modeling  formulation. 
The  case  studies  are  presented  in  Section  4.  Finally,  the  conclusions 
are  provided  in  Section  5. 
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2.  Large-scale  photovoltaic  systems 

The  electrical  structure  of  a  large  scale  PV  system  is  shown  in 
Fig.  1.  The  PV  power  generation  model  should  be  drastically  derived 
as  a  function  of  solar  irradiation  and  surface  temperature  of  solar 
panels  in  order  to  determine  the  solar  power  generation  more 
efficiently  [24  .  The  equivalent  circuit  of  the  PV  module  is  shown  in 
Fig.  2. 

The  I—V  characteristic  of  a  PV  cell  is  expressed  as: 

f  !  I  fr»rn(VpV +  Rs  X  IpV)  ll  VpV  +  X  IPV 

fpv  =  /  ph  -  Is  [exP(^ - Vtxa - j  -  lj - ^ -  (1) 


Ipv 


where 


Vt 


NsxkxTj 

9 


(2) 


Both  the  PV  output  current  and  saturation  current  change  with 
the  solar  radiation  and  panel  temperature.  Thus,  the  PV  internal 
equivalent  current  source  can  be  expressed  as: 


/Ph  =  (/sc  +  K,  x  AT)  x  ^  (3) 

The  reverse  saturation  current  can  be  also  represented  as 
follows: 


Is  —  hr 


exp 


9  X  Eg 
a  x  k 


(4) 


The  maximum  output  power  of  a  PV  cell  is  then  expressed  by 


P m.ix  —  Vp\!  X  Ip 


,rPvx(lpt-i,[ex  P(t'p''^:Wl-1 
VpV  +  Rs  x  Ip 

Rsh 


(5) 


The  modeling  of  a  solar  module  is  discussed  in  more  details  in 
Ref.  [33],  Seven  large-scale  PV  systems,  with  nominal  power  output 
from  100  kW  to  2500  kW,  are  used  in  order  to  calculate  the  overall 
system  reliability.  The  PV  module  and  inverter  characteristics 
which  listed  in  Table  1,  are  the  same  in  all  systems. 


Table  1 

The  PV  module  and  inverter  characteristics. 


Inverter  (100  kW) 

PV  module  (230  W) 

Vmpp.min  =  450  V 

/sc  =  8.24  A 

Vmpp.max  =  820  V 

Vmpp  =  30.2  V 

Vmax  =  1000  V 

fi,  =  3.3  mA  -C-1 

/DC,max  =  235  A 

Voc  =  37.2  V 

/mpp  =  7.60  A 

Hv  =  -120  mV  °C 1 

To  ensure  that  no  inverse  current  is  running  in  a  string,  the  PV 
module  strings  are  connected  to  the  inverter  by  means  of  a  pro¬ 
tection.  The  protection  can  be  a  blocking  diode  (as  shown  in  Fig.  1 ), 
a  fuse  or  a  circuit  breaker.  After  the  string  protection,  the  DC  switch 
helps  to  disconnect  the  PV  field  for  maintenance,  even  under  solar 
irradiation.  After  the  DC  switch,  the  charge  controller  regulates  the 
power  coming  from  the  solar  panels  to  the  batteries  and  prevents 
the  batteries  from  over-charging.  Then,  the  battery  system  stores 
the  generated  energy  for  emergency  use  at  night  and  on  cloudy 
days.  To  convert  from  DC  to  AC,  inverters  are  added,  which  are 
protected  from  lightning  by  surge  protection  devices  (SPDs).  To 
protect  the  AC  lines,  as  per  normal  electrical  design  practice,  an  AC 
circuit  breaker  is  added  after  the  inverter.  Ultimately,  the  PV  system 
is  connected  to  the  transformer,  though  that  was  not  done  in  this 
study.  Table  2  reports  the  number  of  components  in  each  PV  sys¬ 
tem,  which  increased  as  the  system’s  nominal  power  output 
increased.  Note  that  in  real-world  installations,  some  devices  could 
be  eliminated.  In  order  to  focus  only  on  failures  of  electrical  and 
electronics  components,  the  failure  of  cables  was  not  considered. 
Therefore,  the  analysis  assumed  a  flawless  system  installation,  as 


PV  Module  s,rin8 
Protection 


System 


Differential  , 

AC  Circuit  [ 

Switch  Breaker  j  Transformer 


Fig.  1.  Electrical  structure  of  the  large  scale  PV  system. 
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Table  2 

Number  of  components  per  each  PV  system. 


Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

PV  modules 

437 

874 

2166 

4351 

6517 

8702 

10,868 

String  protection 

23 

46 

114 

229 

343 

458 

572 

DC  switch 

3 

6 

15 

27 

42 

57 

72 

Inverter 

1 

2 

5 

9 

14 

19 

24 

AC  circuit  breaker 

1 

2 

5 

9 

14 

19 

24 

Grid  protection 

AC  switch 

Differential  circuit  breaker 
Connector  (couple) 

1 

1 

1 

874 

1 

1 

1 

1748 

1 

1 

1 

4332 

1 

1 

1 

8702 

1 

1 

1 

13,034 

1 

1 

1 

17,404 

1 

1 

1 

21,736 

Battery  system 

Charge  controller 

16 

1 

30 

1 

76 

1 

150 

1 

224 

1 

298 

1 

372 

1 

well  as  flawless  SPDs  with  no  failure  rates.  The  reliability  of  all 
seven  PV  systems  was  analyzed  over  one  year  of  operation  and  over 
20  years  operation  with  8.5  h  of  operation  per  day.  Therefore,  the 
failure  rates  for  all  components  are  failures  hour-1. 

The  batteries  used  in  this  study  came  from  the  Rolls— Surrette 
factory.  All  PV  systems  used  an  identical  battery  model,  Ah  and 
voltage,  12  CS  IIP,  475  Ah  and  12  V,  respectively.  Many  installation 
methods  can  calculate  the  number  and  structure  of  batteries  in  PV 
systems  [34];  however,  in  this  study,  it  was  assumed  that  each 
battery  system  had  two  days  of  reserve  capacity.  It  is  also  important 
to  point  out  that  considering  the  overall  PV  power  output  as  a 
reserve  capacity  is  not  economical.  Therefore,  the  analysis  consid¬ 
ered  only  the  emergency  reserve  needed  in  PV  systems,  which  was 
assumed  to  be  5%  of  the  PV  power  output  for  all  PV  systems.  The 
battery  system  consisted  of  two  batteries  in  series,  each  with  a 
nominal  voltage  of  12,  for  a  total  required  voltage  of  24  and  M 
batteries  in  parallel. 

For  instance,  the  method  used  to  determine  the  number  of 
batteries  in  the  battery  bank  for  a  200  kW  PV  system  is  described  as 
follows: 

1.  Determining  the  emergency  reserve  needed 
(5%  x  200  x  103  W  =  10,000  W). 

2.  Defining  the  number  of  reserved  days  (2  days). 

3.  The  third  step  is  to  determine  the  watt  hours  needed  during  periods 
of  little  or  no  sunshine  (10,000  W  x  2  days  x  8.5  h  =  170,000  W  h). 

4.  Determining  the  amp  hours  needed  in  the  periods  of  no  sun¬ 
shine  from  dividing  the  watt  hour  in  step  three  by  the  battery 
system  voltage  (170,000  W  h/12  V  =  14  250  Ah). 

5.  Dividing  the  amp  hour  calculated  in  the  step  four  by  the  battery 
system  amp  hour  in  order  to  determine  the  total  number  of 
batteries  (14,250  Ah/475  Ah  =  30).  Then,  divide  this  by  2 
in  order  to  determine  the  total  number  of  batteries  in  parallel 
(30/2  =  15). 

This  steps  can  be  summarized  as  follows: 

(200  kW  x  2  days  x  8.5  h  x  5%)/12  V 
N  '  1 - 475Ah - L - ’  30  (6) 

M  -  30/2  =  15  (7) 

3.  Basic  reliability  concepts  and  mathematics 

3.3.  General  reliability  function 

Reliability  is  defined  in  quantitative  terms  as  the  probability  of  a 
component  performing  its  function  adequately.  The  reliability 
function  is  the  probability  that  a  system  will  be  successfully  oper¬ 
ating  without  any  failure  in  the  given  time  t: 


R(t)  =  P(T  >  t)  (8) 

The  unreliability  or  failure  probability  can  be  also  expressed  as: 

F(t)  =  P(T  <  t)  =  1  -  R(t)  (9) 

Equations  (8)  and  (9)  with  a  density  function  f{t)  could  be 
written  as: 

R(t)  =  j  f(t)dt  (10) 


F(t)  =  j  f(t)dt  (11) 

The  mean  time  to  failure  (MTTF)  of  a  component  can  be 
expressed  as: 

MTTF  =  J  t  xf(t)dt  =  J  R(t)dt  (12) 

o  o 

Equation  (12)  represents  the  reliability  of  a  component.  In  the 
real  world,  systems  are  complex  and  consist  of  a  large  number  of 
components  that  may  be  connected  in  series,  in  parallel  or  in  a 
combination  of  series  and  parallel.  When  connected  in  series,  the 
failure  of  one  component  interrupts  the  overall  system,  whereas  all 
components  must  fail  in  the  parallel  system  to  interrupt  the  overall 
system. 

The  reliability  of  a  series  system  with  n  components  can  be 
calculated  as: 

Rsystem  =  P[E  1  H  3y  Pi  ■  ■  ■  O  F,-J  (13) 

If  the  n  components  are  independent,  thus: 

Rsystem  =  P(Ei)  X  P(E2)-  X  P(E;)  (14) 

Assuming  the  system  is  not  repairable  [32],  the  overall  system 
reliability  can  always  be  derived  by  Boolean  techniques.  Thus, 
system  reliability  performance  can  be  expressed  as  a  function  of 
components  reliability  as  follows: 

^system  =  R\  X  R2  X  ...Rn  (15) 

Thus 

n 

Rsystem  =  J  j  Rj  (16) 

i  =  1 

where  R;  =  P(E >■)  is  the  reliability  of  component  i. 

The  system  reliability  of  x  series  units  with  M  parallel  compo¬ 
nents  in  each  unit  can  be  obtained  using: 

Rsystem  =  1  -  (1  -  RX)M  (17) 

A  detailed  discussion  of  the  existing  technique  can  be  found  in 
Refs.  [38-40], 

The  reliability  data  were  taken  from  Refs.  [18]  and  [21—37],  The 
failure  rates  of  the  battery  system  and  the  charge  controller  were 
0.08  and  0.04  for  2  years,  respectively  [36];  hence,  the  failure  rates 
were  calculated  as: 
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^Battery  =  0.08/2/365/8.5  =  12.89  X  10-6  failures  fT1  (18) 

^charge-controller  =  0.04/2/365/8.5  =  6.44  x  ICG6  failures  h_1 

(19) 

All  failure  rates  appear  in  Table  3. 

3.2.  Fault  tree  method 

The  fault  tree  method  (FTA)  is  the  most  useful  tool  in  analyzing 
the  risk  and  reliability  of  causal  systems.  FTA  is  a  graphical  design 
method  in  which  failures  are  defined  more  easily  than  non¬ 
failures.  The  focus  is  usually  on  a  failure  appearing  at  the  top  of 
the  fault  tree  diagram.  One  of  the  major  advantages  of  FTA  is  that 
it  can  predict  the  most  system  failures  in  a  system  breakdown.  FTA 
attempts  to  associate  failure  processes  of  logic  diagrams,  which 
show  the  state  and  behavior  of  the  system.  The  top  event  defines 
the  failure  mode  of  the  system  or  its  function,  which  is  then 
analyzed  in  terms  of  failure  modes  of  its  components  and  influ¬ 
ence  factors  [41],  The  top  event  of  a  fault  tree  should  be  consid¬ 
ered  carefully.  If  it  is  selected,  the  analysis  generally  becomes 
unmanageable,  and  its  success  is  not  guaranteed.  To  clarify  the 
advantages  of  proposed  approach,  a  simplified  fault  tree  of  a 
distribution  system  with  a  probability  of  Pa,  a  protection  system 
with  a  probability  of  PB,  restoring  upstream  customers  with  a 
probability  of  Pc,  and  restoring  downstream  customers  with  a 
probability  of  PD,  is  given  in  Fig.  3.  Each  combination  of  successful 
and  unsuccessful  system  responses  corresponds  to  a  path  on  the 
FTA,  and  the  probability  of  any  outcomes  are  calculated  by 
multiplying  the  associated  probabilities  of  each  individual  system 
response  [42], 

A  fault  tree  also  is  assigned  by  "gates”  that  present  the  re¬ 
lationships  between  their  input  and  output  events.  A  detailed  dis¬ 
cussion  of  fault  tree  techniques  is  outside  the  scope  of  this  study 
but  can  be  found  in  Ref.  [43], 

Fig.  4  shows  a  fault  tree  for  the  PV  system  illustrated  in  Fig.  1. 

A  system  may  have  more  than  one  top  event,  as  shown  in  Fig.  4. 
The  top  event  appears  in  a  box  that  represents  the  failure  event 
under  investigation.  For  instance,  Energy  reduced  if  PV  irradiated 
was  determined  as  a  top  event  in  this  study. 

The  following  basic  symbols  are  used  in  Fig.  4  to  represent  the 
relationship  between  the  top  event  and  lower  level  events: 


Fault  Protection  Upstream  Downstream 

occurs  system  restoration  restoration 


Probability 


Succeeds 


Pd 

Pd 

Pd 

Pd 

Pd 

Pd 

Pn 


PaPbPcPd 
PaPbPcP  d 
Pa?bP  cPd 
Pa?bP  c?  d 
PaP  bPcPd 
PaP  bPcP  d 
P  aP  bP  cP  d 
PaP  bP  cP  d 


Fig.  3.  A  simplified  fault  tree. 


•  Circle:  A  circle  is  a  base  event  in  the  tree  and  does  not  require 
any  gates  or  events  as  inputs. 


3.3.  Minimal  cut  sets  evaluation  using  FTA 

In  this  section,  a  mathematical  model  for  cut  sets  is  presented 
using  FTA.  Such  methods  significantly  help  to  convert  fault  trees 
into  Boolean  models  and  mathematical  equations.  By  definition,  a 
minimal  cut  set  causes  the  system  to  be  unavailable  because  of 
component  failures.  If  all  components  are  unavailable,  a  minimal 
cut  set  will  cause  the  top  event  to  occur.  If  all  components  are 
available,  then  the  top  event  will  not  occur.  The  minimal  cut  set 
for  any  fault  tree  is  finite  and  can  be  achieved  easily.  For  instance, 
the  three  component  minimal  cut  sets  show  that  all  three  com¬ 
ponents  must  fail  in  order  for  the  top  event  to  occur.  For  an  n 
component  minimal  cut  set,  all  n  components  must  be  unavailable 
[43]. 

The  minimal  cut  set  expression  for  the  top  event  can  be  obtained 
from  Fig.  4  as: 


•  OR  gate:  The  event  above  the  gate  will  occur  if  at  least  one 
combination  of  the  input  events  (the  events  shown  below  the 
gate)  exist. 

•  Rectangle:  The  rectangle,  the  main  block  of  FTA,  represents  a 
negative  event.  It  is  located  at  the  top  of  the  tree  or  throughout 
the  tree  to  denote  other  events  of  interest. 


Table  3 

Component  failure  rates. 


Component 

Failure  rate  (10  6  failures  h  ’) 

Reference 

PV  modules 

0.0152 

[181 

String  protection 

0.313 

[35J  Sect.  6-2 

DC  switch 

0.2 

[35J  Sect.  22-1 

Inverter 

40.29 

[21) 

AC  circuit  breaker 

5.712 

[35]  Sect.  14-5 

Grid  protection 

5.712 

[35]  Sect.  14-5 

AC  switch 

0.034 

[35]  Sect.  14-1 

Differential  circuit  breaker 

5.712 

[35]  Sect.  14-5 

Connector  (couple) 

0.00024 

[35]  Sect.  17-1 

Battery  system 

12.89 

[36[ 

Charge  controller 

6.44 

[36[ 

Top  event  =  (((((((((PV  +  CON  +  BD)  +  DCS)  +  CC) 

+  (SPD  +  BS))  +  INV)  +  (CBac  +  SPD))  +  GP) 

+  ACS)  +  DCB) 

(20) 

It  is  also  assumed  that  SPDs  are  flawless  and  have  no  failure 
rates;  therefore,  Equation  (20)  is  reduced  to: 


Top  event  =  PV  +  CON  +  BD  +  DCS  +  CC  +  BS  +  INV  +  CBac 
+  GP  +  ACS  +  DCB 


(21) 


Equation  (21)  indicates  that  the  fault  tree  in  this  study  had  11 
minimal  cut  sets.  The  analysis  of  the  FTA  can  be  categorized  into  the 
two  kinds  of  qualitative  and  quantitative  analysis.  The  qualitative 
analysis  is  obtained  from  minimal  cut  sets  (the  combinations  of 
events  which  can  cause  system  failure)  and  quantitative  fault  tree 
evaluation  deals  with  calculating  the  top  event  probability  ac¬ 
cording  to  the  bottom  events  [44],  The  probability  of  failure  for  top 
event  can  be  calculated  on  qualitative  analysis  as  follows  [41]: 
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Energy  null  or  reduced 
if  PV  irradiated(Top  Event) 


n 

^TOP  =  X!^Mcsi  “  y^MCSinMCSi  +  ^MCSinMCSjnMCSt - 

i  =  1  i<j  i<j<k 

1  m 

+  (-l )  p  n  MCS, 

i  =  l 

The  probability  of  occurrence  of  MCS,  can  be  expressed  as:  (22) 

PmCS,  =  Pb  1  x  Pb2\PB1  x  OPB2  X 

x  ^BnlPfil  OPB2  O  ■  PlPBn-l  (23) 

It  was  assumed  that  the  components  are  independent.  Thus, 

Pmcs,.  =  n  PBi  (24) 

i  =  1 

Using  probability  theory  and  the  minimal  cut  set,  the  probability 
of  failure  of  the  top  event  can  be  obtained  from  Equation  (21),  and 
the  reliability  probability  can  be  expressed  as: 


Pr(Top  event)  =  Pr(E^  +  E2...En)  (25) 

1  - Pr(Top event)  =  [1  - Pr(E 1 )]  x  [1  -Pr(l  - E2)]. . . [1  - Pr(l  - £„)] 

(26) 

It  is  assumed  that  the  components  are  independent,  so  the  total 
system  reliability  can  be  calculated  as: 

Rial  =  ]  [ 
i  =  i 

=  R(PV)  x  R(CON)  x  R(BD)  x  R(DCS)  x  R( CC)  x  R(BS) 
x  R(INV)  x  R(CBac)  x  R(GP)  x  R( ACS)  x  R(DCB) 

(27) 

Using  exponential  distribution  and  Equation  (16),  the  total 
system  reliability  can  be  shown  as: 
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Table  4 

Total  component  reliability  for  the  PV  systems  for  a  period  of  one  year  of  operations  [in  %]. 


Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

PV  modules 

97.9603 

95.9622 

90.2899 

81.4497 

73.5409 

66.3405 

59.8988 

String  protection 

97.7913 

95.6313 

89.5204 

80.0613 

71.6712 

64.0981 

57.3809 

DC  switch 

99.8140 

99.6284 

99.0736 

98.3386 

97.4276 

96.5250 

95.6307 

Inverter 

88.2497 

77.8801 

53.5262 

32.4653 

17.3775 

9.3015 

4.9787 

AC  circuit  breaker 

98.2435 

96.5178 

91.5205 

85.2576 

78.0281 

71.4117 

65.3564 

Grid  protection 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

AC  switch 

99.9895 

99.9895 

99.9895 

99.9895 

99.9895 

99.9895 

99.9895 

Differential  circuit  breaker 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

98.2435 

Connector  (couple) 

99.9349 

99.8699 

99.6780 

99.3541 

99.0342 

98.7125 

98.3946 

Battery  system 

99.9996 

99.9992 

99.9981 

99.9962 

99.9942 

99.9923 

99.9904 

Charge  controller 

98.0218 

98.0218 

98.0218 

98.0218 

98.0218 

98.0218 

98.0218 

Table  5 

Total  component  reliability  for  the  PV  systems  for  a  period  of  20  years  of  operations  [in  %].  Note  that  0%  does  not  mean  that  the  overall  PV  system  is  failed,  it  means  that  at  least 
one  component  of  the  PV  system  is  failed. 


Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

PV  modules 

66.2218 

43.8532 

12.9654 

1.6512 

0.2141 

0.0273 

0.0035 

String  protection 

63.9737 

40.9264 

10.9255 

1.1707 

0.1279 

0.0137 

0.0015 

DC  switch 

96.3455 

92.8245 

83.0149 

71.5288 

59.3796 

49.2939 

40.9213 

Inverter 

8.2085 

0.6738 

0.0004 

0.0000 

0.0000 

0.0000 

0.0000 

AC  circuit  breaker 

70.1574 

49.2205 

16.9968 

4.1177 

0.6999 

0.1190 

0.0202 

Grid  protection 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

AC  switch 

99.7893 

99.7893 

99.7893 

99.7893 

99.7893 

99.7893 

99.7893 

Differential  circuit  breaker 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

70.1574 

Connector  (couple) 

98.7069 

97.4305 

93.7525 

87.8455 

82.3574 

77.1684 

72.3473 

Battery  system 

99.9923 

99.9846 

99.9616 

99.9232 

99.8848 

99.8464 

99.8080 

Charge  controller 

67.0587 

67.0587 

67.0587 

67.0587 

67.0587 

67.0587 

67.0587 

Rtot  =  exp  ^  mfkitij  (28) 

where  m,-  is  the  total  number  of  components  (i.e.  8702  PV  modules 
for  2  MW  PV  system),  A,-  is  the  failure  rate  of  component  i,  n  is  the 
total  number  of  different  components,  and  t  is  the  study  time  of  the 
reliability  analysis.  For  instance,  the  total  reliability  of  the  PV 
modules  for  a  2  MW  system  (8702  PV  modules)  for  one  year 
(365  days  x  8.5  h  =  3102.5  h)  is: 

Rpv.tot  =  exp(-  1  x  8702  x  0.0152  x  10  6  x  3102.5)  =  0.663405 

(29) 


declined  quickly.  For  a  200  kW  system,  the  PV  modules  were 
43.8532%  reliable,  while  the  inverter  was  only  0.6738%  reliable.  For 
a  2.5  MW  system,  the  PV  modules  had  only  a  0.0035%  probability  of 
operating  correctly,  while  the  inverter  was  not  reliable,  with  a  0% 
probability  of  operating  without  failures  (see  Table  5). 

A  reliability  of  0%  means  that  there  is  at  least  one  component 
with  a  failure,  which  induces  the  failure  of  the  overall  PV  system. 

The  overall  system  reliability  for  one  year  and  for  20  years  of 
operation  is  reported  in  Table  6. 

To  obtain  useful  information  from  this  study,  the  Fussel— Vesely 
method  was  used.  This  method  accounts  for  the  effect  of  each 
single  component  on  the  overall  system  reliability,  which  can  be 
calculated  as: 


The  total  component  reliability  for  the  PV  system  over  one  year 
and  over  20  years  of  operation  was  calculated  using  Equation  (28), 
(see  Tables  4  and  5).  Additionally,  the  failure  rate  for  the  battery 
system  that  consists  of  series  and  parallel  paths  was  determined 
using  Equation  (17);  for  instance,  the  total  reliability  of  the  battery 
system  for  a  100  kW  PV  system  is  represented  as: 


(31) 


For  instance,  the  Fussel— Vesely  of  the  string  protection  for  a 
100  kW  system  (23  string  protections)  for  20  years  of  operation 
(365  days  x  20  years  x  8.5  h  =  62,050  h)  is: 


^Battery-system  =  1  -  (l  -  R2)&  =  1  328  X  10  9  (30) 

4.  Simulation  and  results 

By  substituting  the  failure  rates  listed  in  Table  3  into  Equation 
(28),  the  total  component  reliability  of  PV  systems  for  one  year  and 
20  years  of  operation  was  estimated  (See  Tables  4  and  5).  The 
component  reliabilities  decreased  as  the  PV  power  output 
increased;  for  instance,  after  one  year  for  a  200  kW  system,  the  PV 
modules  had  a  95.9622%  probability  of  operating  without  failures, 
while  the  inverter  had  only  a  77.8801%  probability.  For  a  2.5  MW 
system,  the  PV  modules  had  a  59.8988%  probability  of  operating 
without  failures,  while  the  inverter  had  only  a  4.9787%  probability 
(see  Table  4).  However,  for  20  years  of  operation,  the  reliability 


FV  =  [l  -  (exp(-  lx  23  x  0.313  xl0~6x  62050))]  =0.36023 

(32) 

The  Fussel— Vesely  method  reveals  the  impact  of  the  PV  mod¬ 
ules,  string  protections,  DC  switches,  inverters,  AC  circuit  breakers, 

Table  6 

The  overall  system  reliability  for  a  period  of  one  and  20  years  of  operations  [in  %]. 
Note  that  0%  does  not  mean  that  the  overall  PV  system  is  failed,  it  means  that  at  least 
one  component  of  the  PV  system  is  failed. 


Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

Reliability 
(in%,  1  year) 

78.3716 

64.9282 

36.9896 

16.6818 

6.5229 

2.5457 

0.9954 

Reliability 

(in%,  20  years) 

0.7641 

0.0177 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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Table  7 

Results  of  Fussel— Vesely  for  one  year  of  operations. 


Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

PV  modules 

2.0397 

4.0378 

9.7101 

18.5503 

26.4591 

33.6595 

40.1012 

String  protection 

2.2087 

4.3687 

10.4796 

19.9387 

28.3288 

35.9019 

42.6191 

DC  switch 

0.1860 

0.3716 

0.9264 

1.6614 

2.5724 

3.4750 

4.3693 

Inverter 

11.7503 

22.1199 

46.4738 

67.5347 

82.6225 

90.6985 

95.0213 

AC  circuit  breaker 

1.7565 

3.4822 

8.4795 

14.7424 

21.9719 

28.5883 

34.6436 

Grid  protection 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

AC  switch 

0.0105 

0.0105 

0.0105 

0.0105 

0.0105 

0.0105 

0.0105 

Differential  circuit  breaker 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

1.7565 

Connector  (couple) 

0.0651 

0.1301 

0.3220 

0.6459 

0.9658 

1.2875 

1.6054 

Battery  system 

0.0004 

0.0008 

0.0019 

0.0038 

0.0058 

0.0077 

0.0096 

Charge  controller 

1.9782 

1.9782 

1.9782 

1.9782 

1.9782 

1.9782 

1.9782 

Table  8 

Results  of  Fussel— Vesely  for  20  years  of  operations. 

Power  (kW) 

100 

200 

500 

1000 

1500 

2000 

2500 

PV  modules 

33.7782 

56.1468 

87.0346 

98.3488 

99.7859 

99.9727 

99.9965 

String  protection 

36.0263 

59.0736 

89.0745 

98.8293 

99.8721 

99.9863 

99.9985 

DC  switch 

3.6545 

7.1755 

16.9851 

28.4712 

40.6204 

50.7061 

59.0787 

Inverter 

91.7915 

99.3262 

99.9996 

100 

100 

100 

100 

AC  circuit  breaker 

29.8426 

50.7795 

83.0032 

95.8823 

99.3001 

99.8810 

99.9798 

Grid  protection 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

AC  switch 

0.2107 

0.2107 

0.2107 

0.2107 

0.2107 

0.2107 

0.2107 

Differential  circuit  breaker 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

29.8426 

Connector  (couple) 

1.2931 

2.5695 

6.2475 

12.1545 

17.6426 

22.8316 

27.6527 

Battery  system 

0.0077 

0.0154 

0.0384 

0.0768 

0.1152 

0.1536 

0.1920 

Charge  controller 

32.9413 

32.9413 

32.9413 

32.9413 

32.9413 

32.9413 

32.9413 

Total  component  reliability  for  PV  systems  after  one  year  of  operations 
100 

80 
60 
40 
20 
0 

FusseleVesely  relative  importance  measures  for  a  period  of  one  years  of  operations 


Reliabilities  after  one  years  of  operations 


100 


1  2  3  4  5  6  7 


Total  component  reliability  for  PV  systems  after  twenty  years  of  operations 
100 

80 
60 
40 
20 
0 

FusseleVesely  relative  importance  measures  for  a  period  of  twenty  years  of  operations 
100 
80 
60 
40 
20 
0 


1  2  3  4  5  6  7 

Reliabilities  after  twenty  years  of  operations 


1  2  3  4  5  6  7 


DCb  I  I  CON  battery  charge-controller 


Fig.  5.  The  all  results  of  reliability  for  seven  PV  systems. 
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Table  9 

Critical  component  priorities. 


Priority 

Component 

1 

Inverter 

2 

String  protection 

3 

PV  modules 

4 

AC  circuit  breaker 

5 

DC  switch 

6 

Charge  controller 

7 

Grid  protection 

7 

Differential  circuit  breaker 

9 

Connector  (couple) 

10 

AC  switch 

11 

Battery  system 

grid  protection,  AC  switches,  differential  circuit  breaker,  connec¬ 
tors,  battery  systems  and  charge  controllers  on  the  overall  system 
reliability  of  large-scale  PV  systems. 

The  Fussel— Vesely  results  for  one  year  and  for  20  years  of 
operation  appear  in  Tables  7  and  8. 

Fig.  5  displays  all  of  the  results  of  this  study  for  seven  PV  system 
tests,  numbered  1  to  7. 

5.  Conclusion 

Photovoltaic  systems  are  fairly  reliable  and  usually  can  be 
trusted  to  operate  without  any  failures.  However,  like  any  other 
system,  they  may  fail,  so  the  effects  of  failures  should  be  calculated. 
This  paper  described  a  method  for  evaluating  the  reliability  of  the 
components  in  large-scale  PV  systems.  The  ranking  of  the  most 
critical  components  appears  in  Table  9.  The  component  that  fails 
the  most  must  be  able  to  be  replaced  or  repaired  quickly  after  a 
failure.  PV  systems  could  have  a  high  rate  of  availability  under 
frequent  maintenance,  but  frequent  maintenance  is  not  an  optimal 
solution.  Maintenance  strategies  can  be  optimized  to  reduce  the 
associated  costs.  It  is  necessary  to  point  out  that  although  the  scope 
of  this  paper  was  based  on  unrepairable  components,  the  proposed 
method  can  be  further  considered  with  the  repair  time  in  order  to 
evaluate  the  system  economic  losses.  In  addition,  Future  work  will 
focus  on  which  types  of  maintenance  are  appropriate  for  the 
components  of  large-scale  PV  systems. 
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